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ABSTRACT 

The hierarchical growth of dark matter haloes, in which galaxies are hosted, has been 
studied and modeled using various approaches. In this paper we use a modified version 
the Sheth & Lemson algorithm for a A cold dark matter power spectrum, and model 
the growth of a Milky- Way sized halo with microsolar mass resolution, corresponding 
to the typical Jeans mass for a dark matter Weakly Interacting Massive Particle with 
mass of 100 GcV. We then compute the unevolved subhalo mass function and build-up 
a Milky- Way halo placing and evolving its satellites. This subhalo population is used 
to study the 7-ray emission from dark matter annihilation. In this case, the subhaloes 
which populate the host halo have been computed considering only progenitor haloes 
accreted by the main branch of the tree, so as to correctly treat the embedding of 
sub-subhalocs inside subhaloes. Each subhalo will indeed host at the present-time 
sub-subhaloes accreted when it was an isolated system. In order to compute the sub- 
subhalo population of a Milky- Way dwarf galaxy, like Draco, and to study its 7-ray 
emission, we first estimate the Draco virial mass at merging redshift z m and then we 
run the merger tree from z m following the halo down to the dark matter Jeans mass. 
We then study the effect on the Fermi-LAT (GLAST) detectability for both subhaloes 
in the Milky- Way and in Draco, and we show how subhaloes cannot be responsible 
for the boost factor needed for detection. 
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1 INTRODUCTION 

In the standard scenario of structure formation galaxies re- 
side in massive Dark Matter (DM) haloes, where baryons 
can shock, cool and eventually form stars. Haloes form by 
gravitational instability, starting from some density fluctu- 
ation field S(x). Specifically, a dark matter halo with mass 
M oc R 3 , forms when the linear density fi eld, smoothed on 
scale R, exceeds some critical threshold S c ( Bond et al.lll99ll : 
ICarroll et al1ll992l : lLacev fc Coielll993l : lEke et alJll99a\ 

In a Cold Dark Matter (CDM) framework, the galaxy 
formation process is hierarchical along cosmic time: small 
systems collapse earlier, in a denser universe, and later 
merge to form larger haloes. In particular, haloes of- 
ten grow as a consequence of repeated merging events 
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with smaller satellites. Present-day surviving satellites form 
the so-called subh alo (or substructure) population of a 
given host system (IMoore et al.lll999l: ISpringel et al.ll200ll: 



Gap et al.l |2004|; Ivan den Bosch et al.l 120051: iGiocoli et all 



2008al : IZentner et alJl2005F 



Dark matter haloes and subhaloes pro vide the envi- 
ronment in which galaxies form and evolve (IBullock etal 



200d : ISomervillell2002l : iKravtsov et~alll2004l ; IVale fc Ostriker 



200a ). Therefore, understanding the assembly histories of 
dark matter haloes is the first step towards the compre- 
hension of the more complex processes involved in galaxy 
formation. 

In order to study the framework of structure formation, 
two different approaches have co mmonly been used. T h e first 
is to run TV-body sim ulations (jSpringel et all l200ll . 120051 : 
Dic mand et al ] |2007al lbh. These are very powerful tools, that 
can be used to reproduce the collapse of dark matter haloes 
with high mass and force resolution, both on galaxy and 
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galaxy-cluster scales. They allow to follow gravitational col- 
lapse up to its fully non-linear evolution. On the other hand, 
simulations are computationally expensive, and nonetheless 
cannot cover the full spectrum of masses relevant to struc- 
ture formation. The second approach is to use some analyt- 
ical modeling, which allows a detailed study of the merging 
history of haloes over an arbitrarily large mass range, un- 
der suitable simplification of the problem. This is the case of 
the Press fc Schechter formalism llLacey fc Coldl 19931 : IShethl 
1 19981 ; ISheth fc Lemsonll 19991 ; IShethll2003l )~ 

As underlined before, haloes collapse on a certain 
scale once the linear density contrast smoothed on that 
scale exceeds some threshold value. The nonlinearities in- 
troduced by these virialized objects do not affect the col- 
lapse of overdense regions on larger scales. This simple as- 
sumption leads to the derivation of the global mass func- 
tion of dark matter haloes that is in fair agreemen t with 
that found in iV-body simulations (|Lacev fc Cold 1 19931 : 
ISomerville et all |2000| ; ISheth fc Tormenlll999l . |2002|) . Dur- 
ing the last decades, an extension of thi s theory was made 
by several a uthors (jLacev fc Cold 1 19931 ; iBond et all Il99ll ; 
iBowerl Il99ll ) with the aim of computing the probability 
that a halo of mass m, at redshift 2, belongs to a given 
halo of mass Mp at r edshift 20 < 2 ()Lacev fc Cold 1 19931 : 
ISheth fc Tormenll2002l) . This quantity was named "condi- 
tional" or "progenitor" mass function, f(m, z\Mo, 2o)dm. 
This allowed di fferent authors to estim ate quantities such 
as the merger (jNeistein fc Dekell 120081) and creation rate 
llPertival fc Milled 1 19991 ; IPercival et al.ll2000l : iMoreno et al.1 
120071 ). and the formation time distribution of dark matter 
haloes (|Sheth fc Tormenll2004 iGiocoli et al.ll2007l ). 

The extended Press-Schechter formalism can be numer- 
ically implemented to produce stochastic realizations of the 
merging history tree of haloes of any mass; these Monte 
Carlo merging histories can have, theoretically, arbitrary res- 
olution in mass and time, and their result is the full popu- 
lation of progenitor haloes at all times, for some final halo. 
The only thing they do not provide is spatial information on 
the haloes themselves. Therefore, when the focus is on the 
statistical properties of haloes and not on their position or 
internal structure, Monte Carlo merger trees have the ad- 
vantage - over iV-body simulation - of virtually unlimited 
mass resolution. Implementing a code to reach the required 
resolution is however not straightforward, and caution must 
be used in order to preserve consistently with the theory. 

In this paper we will build a spherical collapse Monte 
Carlo merger tree with arbitrary mass and time-step reso- 
lution, and will consider a ACDM-power spectrum. We will 
populate a Milky-Way sized halo with subhaloes with masses 
as small as 10 _6 Mq. This value correspond to the typical 
Jeans mass for the a CDM Weakly Interacting Massive Par- 
ticle (WIM P) particle with m D M = WOGeV (|Green et al.1 
I2004 I2005I ). Such a value for the minimum mass can ac- 
tually vary between 10~ 12 and 10~ 4 Mp) depend ing on the 
underlying particle physics (|Profumo et aLlfeoOfj ). 

We will take into account the progenitor haloes accreted 
by the main branch of our Milky- Way like halo, in order to 
obtain a snapshot of the spatial distribution of substructures 
today. In this way we will select only the first or- 
der of substructures which populate the halo today 
(those which are orbiting in the host halo potential), 
in order to correctly model the radial dependence 



of subhalo properties. In other words, we will avoid 
the bia s of our previous works (see. e.g. IGiocoli et al.1 
towards the small scale structures. In those 
works, the embedding of sub-subhalos within subha- 
los was not treated correctly. Each structure existing 
inside the main halo was indeed considered as orbit- 
ing inside the potential of the host halo itself. In this 
way, an analytical treatment of the effect of subhalos 
on the expected photon flux from DM annihilation 
was possible. However, it was not possible to take 
into consideration the fact that a number of small 
scale structures are indeed sub-substructures, that 
is to say they are orbiting inside the potential of the 
subhalo they belong to, and not inside the poten- 
tial of the host halo. In the present work we will be 
able to separate and study the effect of introducing 
sub-subhalos. With the same partition code we will thus 
compute the subhalo population of subhaloes (i.e. subhaloes 
within subhaloes) , and we will consider the special case of a 
Draco-like satellite. In the second part of the paper we will 
estimate the 7-ray emission from dark matter annihilation, 
in the subhalo population of the Milky- Way and of a Draco- 
like subhalo, and will discuss the different contributions to 
7-ray emissions due to the smooth and clumpy components 
of the considered systems. 

The paper is organized as follows: in Seed] we describe 
the m erger tree technique developed by l|Sheth fc Lemsonl 
Il999l ) and its generalization to a ACDM power spectrum. 
Sec. [3] is dedicated to the study of the merger tree of a 
Milky- Way sized halo, and to compute its mass accretion 
history along cosmic time and its satellite mass function. 
In Sec. U] we model the hierarchical growth of a sample of 
Draco-like satellites until their merging time with the main 
Milky- Way progenitor halo. In Sec. [5] we study the 7-ray 
emission from dark matter annihilation in subhaloes and 
sub-subhaloes and estimate the possibility to be detected 
with the Fermi-LAT telescope. In Sec(S]we discuss our re- 
sults and conclude. 



2 METHOD AND MERGER TREE 

The simplest algorithm for a merger tree uses a binary split. 
In this scenario each halo is split in two haloes at an ear- 
lier epoch. Each of these is in turn divided in other two 
pieces, and so on until all halo masses fa ll below an arbitrar> 



pieces, and so on until all nalo masses ta ll bel ow an arb itrary 
chosen and desired mass resolu t ion m f ([Lacev fc Cole 19931 ; 
ICole fc Kaiserill988l;IColdll99ll :| Kauffmann fc Whitdll993l ). 
However, Somerville fc Kolattl ( 19991 ) showed the failure of 
using a binary merger tree. In this way the first halo which 
is chosen from the Press & Schetchter distribution follows 
the correct probability, but the second one is chosen only in 
order to conserve mass and does not follow the theoretical 
model as expected. This leads to a conditional mass function 
and formation redshift distribution in disagreement with the 
exte nded Press fc Schechter pre dictions. To solve this prob- 
lem, ISomerville fc Kolatd ij 19991 ) proposed a new algorithm 
able to reproduce the conditional mass function quite well, 
at the expenses of using a finely tuned grid of tim e-steps. 

On t he othe r hand , ISheth fc Lemsonl l)l999h used the 
results of IShethl I 1996h and realized that, for white- noise 
initial conditions, mutually disconnected regions are mutu- 
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ally independent. In this case it is possible to split a halo 
into progenitors very efficiently. The so obtained conditional 
mass function is in excellent agreement with the theoretical 
model of the spherical collapse. The great advantage of this 
method is that it is possible to obtain arbitrary high mass 
resolution for any given time-step, generating progenitors in 
a very fast and efficient way. 

In this s ection we will des c ribe the way the method 
developed by ISheth fc Lemsonl (|l999T ) can be generalized 
considering a ACDM power spectrum Using this algorithm, 
we will follow the creation and the assembly history of a 
present-day Milky- Way sized halo. 

Assuming that the dark matter is a WIMP, the Jeans 
mass of the smaller DM halo is given by its free-streaming 
scale ijGreen et al.ll2004l . I2005I ). and for a DM particle mass 
of 100 GeV this smallest halo mass has a typical mass 
mj = 10~ 6 Mq. We therefore extend our merger tree down 
to microsolar mass resolution, in order to study the present- 
day subhalo population down to the dark matter Jeans mass. 

2.1 Poissonian Distribution and Gaussian Initial 
Conditions 

Let us start consid ering an initial Pois sonian distrib ution of 
identical particles. lEpsteinl (|l983l ) and lShethl (|l995l ) showed 
that the probability to find a clump c ontaining N par ticles 
is expressed by the Borel distribution (|Borel E. II 19421 ): 



(Nb)' 



(1) 



where the v ariable s N a nd b are such that N ^ 1 and 
< b < 1. IShethl (119951 ) also showed that this equation 
can be extended to the continuous case, in order to describe 
dark matter halo clustering. In this case the variable 6 has 
a redshift dependence given by: 

6 = 1/(1 + «„(«)) . 

where S c (z) is the critical overdensity threshold predicted 
by the spherical collapse model, decreasing as the universe 
expands. The probability that a randomly chosen particle 
belongs to an A-clump is given by f(N, b) = (l — b)Nr](N, b). 
For large values of N and small values of 5 C , the factorial 
Stirling's approximation gives: 



8 C 



(1 + S C 

5 C 



N 



JV - 1 e -N/(l+S e ) 



(N-iy. 



exp 



2 



(2) 



Since all matter i s in c lump s, it holds the relation 
T,n=i f( N > s <d = l. lShethl (|l995h showed also that Eq. © 
can be obtained starting from an initial Gaussian density 
field with white noise initial conditions (P(k) ~ fc™ with 
n — 0). This equivalence underlines that the Poissonian dis- 
tribution can be thought as the analog ue discrete of the 
white noise Gauss i an po wer spectrum l|Bond et al.l Il99ll : 
lLacev fc Cole]|l993l . [l994l ). 

The critical collapse overdensity decreases with the ex- 
pansion of the universe, so that small systems collapse earlier 
than large ones, i.e., dark matter halo clustering progresses 
hierarchically. An important quantity that describes the hi- 
erarchical growth of the haloes is the conditional distribu- 
tion. It gives the probability that a particle, belonging to a 



clump with N particles at time bo, is part of an n-clump at 
bi < bo, and can be written as: 



/(n,6i|JV,6o) 



bo 



(3) 



where 1 < n < N and < 61/60 < 1 l|Shethlll995h . 

However, for a complete description of the merging his- 
tory tree of an A r -clump at 60, we also need to know the prob- 
ability that, at 61, it is divided in a sample of rij j-clumps 
with k subfamilies (so that ni + ... + tin = k). Recalling 
the conservation of the particle number, X/J=i 3 n i ~ 
the Poissonian Galton- Watson branching process gives the 
following equation for this probability: 

-JV(6i-6 ) 

p(m, rife, 61 1 AT, 60) = 



[Af(6i-6o)f- 



k 

n 

3=1 



v(N,b ) 



(4) 



whe re rj(l,b) is th e Borel distribution with time parameter b 
(see lShethlljl996l ) for more details). Since bo can be related 
to a density, the volume of the A^-clump can be written as: 



Vn,o 



N 



(5) 



n(l + 5 c ,o) ' 

where n denotes the average density of the universe. Eq. Q 
can be thought of as the probability that a region with N 
particles, having density ri(l + 5 c ,o), contains k subregions 
each at average density ri(l + <5 c ,i), where <5 c ,i ^ <5 c ,o- 

The merging history of a present-day halo is de- 
scribed using the E xtended-Press & Schechter formalism 
|Lacev fc Cole! Il993l ) and formulated in terms of its con- 
ditional mass function along consecutive time steps. Let us 
now consider a n n-subclump of the A^c lump, at time 61; in 
Appendix A of ISheth fc Lemsonl (|l999l ) it is shown that, if 
V n ,i = n/ri(l + <5 c ,i) is its associated volume, the remaining 
particles N — n will occupy a volume such that: 

"(i + <y = 77 , (6) 



Vn,0 — Vn, 



where b' 1/(1 + 5' c ) is the unknown quantity to be determined 
and represents the density in the volume Vn-ti- 

The spherical collapse model predicts that a group of 
uncollisionless N particles, with the same mass rriDM, col- 
lapses forming a dark matter M-halo (M = N x 771dm), if 
the smoothed density fluctuation filed on scale R ~ M 1 ^ 3 
exceeds its predicted critical virial value (|Eke et al.lll996l ). 
For each collapsed system we can define its associated mass 



variance as: 



S(M) 



d 3 k(\5{k)\ 2 )W 2 (kR), 



(7) 



where W(kR) and 5(k) represent respectively the Fourier 
transform of the smoothing window function and of the den- 
sity fluctuation field. 

It is possible to show that in the continuous limit 
Eq.s |T} and (g|) give: 



f(y)dp = 2 



du. 



(8) 
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Figure 1. Spherical collapse conditional mass function at six differen t redshi fts, obtained using the SL99 tree. We considered a ACDM 
power spectrum extrapolated down to 10 — 6 Mq as in iGiocoli et al. In each panel the histogram shows the average of 10 4 

realizations of partitioning an Mo-halo. The solid line represents the prediction of the spherical collapse model at the considered redshift. 
The partition was performed both with a single and a multi-step: the distributions are equivalent. 



Taking v = 5 c /vS, this equation describes, at a fixed red- 
shift (and so at a given S c ), the number of collapsed systems 
with mass variance between S and S+dS. The factor 2 takes 
into account the so-called cloud-in- cloud problem, which is 
the possibility that at a given instant some object, which is 
nonlinear on scale M, can be later contained within another 
object, on a larger mass scale. Eq. © describes also the con- 
ditional mass function at time <5 0> i , considering a halo at time 
S Ct o with mass variance So, when v = (5 c ,i — S c .a)/V s — So- 
Before taking into account a ACDM power spectrum, 



let us consider the simpler case of a white-noise power spec- 
trum. We recall that when P(k) ~ fc™, Eq becomes 
M ~ S- (n+3)/3 , that for n = gives M ~ 1/S. Let us 
suppose we want to split an Mo-halo, at time <5 Cj o, in pro- 
genitor haloes at 5 c ,i- To generate the first progenitor we 
draw a random number v\ from the Gaussian distribution, 
Eq and compute its associated mass variance, s\ con- 
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sidering that 



where for this first progenitor m e = Mo (s mE ~ l/m e ) 
and (5 c , me = <5 c ,o- Its physical mass, mi, can be directly 
computed as the inverse of its mass variance. Since for a 
white-noise power spectrum disconnected volumes are mu- 
tually independent, the overdensity of the remaining mass 
m e = Mo — rhi will be given by the volume conservation 
relation in the continuous limit: 

x x ($c,i - 5 cfi ) 

Oc,m E = c .l 7-r- . (10) 

m e /M 

To generate the second halo we draw another number from 
the Gaussian distribution, compute the mass variance from 
the Eq. ©, and hence derive the corresponding mass. The 
remaining mass will now be m e = Mo — rhi — rri2, with the 
corresponding overdensity given always by Eq. (IK)l l. Keep 
going with this procedure it is possible to generate a sample 
of progenitor haloes until the desired m E mass resolution. 

2.2 ACDM Power Spectrum 

In the case of a general power spectrum the algorithm 
should be modified. The assumption that disconnected vol- 
umes are mutually independent does not hold anymore when 
the initial condit i ons d iffer from white-noise. Despite this, 
ISheth fc Lemsonl (|l999h noticed that, when expressed as a 
function of the variance rather than of the mass, all excur- 
sion set quantities are independent of the power spectrum. 
In this framework, each chosen mass m, can be treated not 
as a progenitor having a mass fh, but as a region of volume 
v containing a mass rh, populated by £ objects having all 
the same mass fi, with £ = rh/ [i. 

The number of objects is obtained by requiring that 
they have the same mass variance, that is s(/i) = 1/fh. 
For a scale- free power spectrum P ~ k" (a — n + 3/3) 
£ = m' Q_1 ^ a , and for n = we have f = 1, the region v 
contains exactly one halo, as seen in the previous section. 
For n 7^ and general power spectrum, £ is neither unity 
nor even integer. However, we will show in the next section 
that considering £ = NINT(m//i) (i.e. considering the near- 
est integer to the mass ratio) the progenitor mass functions 
are in excellent agreement with the theoretical prediction at 
all redshifts and down to m e = mj = 

We recall that the ACDM and the white-noise power 
spectrum should satisfy the relation: 

Swn(Afo) = sacdm(Mo) , (11) 

where Mo is the initial mass to be split. This guarantees 
that for the considered initial mass there is one halo in the 
volume Vo for both power spectra. 



3 PARTITION OF A MILKY- WAY SIZE HALO 
WITH MICRO-SOLAR MASS RESOLUTION 

Let us now take into account the case of a ACDM power 
spectrum. The density parameter and mass variance (Qa = 
0.7, Q m = 0.3 and as = 0.772) have been c hosen to agree 
with the recent 3-year WMAP data release (|Spergel et alj 



2007). We have linearly extrapolated the mass variance 
down to the mj = 10 -6 M@ integrating th e power spec- 
trum u sing a top-hat filter in the real space (jGiocoli et al.l 
l2008bl ). In order to have one physical mass both for the 
ACDM and white-noise power spectrum, we should con- 
sider a white-noise power spectrum normalized such that 
Eq. pip holds. For this reason, because we think in term 
of the mass variance than the fisical mass, the mass res- 
olution m t , for a ACDM power spectrum, corresponds to 
the mass resolution ra e , wn for the white-noise one, such that 
Sim(m e , TO ) = s(m € ), with m e , wn > m t . 

3.1 The main branch and the satellite mass 
function 

Let us consider a present-day Milky- Way sized halo (Mo = 
10 12 Mq/H) and suppose we wish to generate a sample of 
progenitors at different redshifts down to m t = mj = 
10" 6 M o , using a ACDM power spectrum. Such mass reso- 
lution corresponds to m £>wn = 2.15 x 10 10 Mq/H. We pro- 
ceed as described in the previous section and obtain the 
conditional mass function plotted in Fig. [T] The number of 
progenitor haloes in the ACDM power spectrum, for a given 
value of the mass variance, was obtained by computing the 
nearest integer of the ratio rh/ [i, where s^ifh) = sacdm(m)- 
In each panel the dotted histogram shows the result of av- 
eraging 10 4 realizations, and the solid line shows the theo- 
retical prediction from the spherical collapse model. At each 
redshift we generated progenitors both with a single and a 
multi-steps technique starting from z = 0. Since the SL99 
method is independent on the time-step, the results from 
the two different ways of progress are in perfect agreement. 

In order to generate the merging history tree along the 
halo main branch, we ran a sample of 10 4 white-noise tree 
realizations for an Mo = 10 12 Mq/Ji initial halo and a mass 
resolution of 0.0215 x Mo. As stressed by SL99 (and demon- 
strated in Fig[T}, this method is time-step independent: to 
obtain a fine enough description of the merging of satellite 
haloes we chose a time resolution AS C = 0.01. We followed 
the main branch (the most massive halo among the progen- 
itors) along consecutive time-steps, storing all the informa- 
tion about t he accreted p r ogenito r haloes (termed "satel- 
lites", as in iGiocoli et al.l (|2008a) ) . until the mass of the 
main branch becomes as small as the mass resolution. We 
call merging redshift z m the most recent redshift before a 
satellite is incorporated in the main halo. 

Finally, the poissonian trees is converted in the ACDM 
one as explained in the previous section. We recall that each 
progenitor halo in the first tree is converted into £ haloes in 
the second, in order to conserve the mass variance in both 
power spectra. 

3.2 Unevolved subhalo population 

IGiocoli et af] (|2008ri) studied the substructure population 
of a present-day Milky- Way halo considering its progenitor 
mass functions at any redshift z > 0, and assuming that all 
progenitor haloes survived and retained their original virial 
mass until redshift zero. This analytical distribution for the 
progenitor mass function was extrapolated down to microso- 
lar mass resolution in order to study the 7-ray emission from 
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Figure 2. Satellite mass function for a present-day Milky- Way 
size halo. The histogram shows the result of 10 4 realizations of 
the merger history tree. The mass of each satellite is rescaled in 
units of the final host halo mass. The solid line shows the slope 
of the least squares fit to the histogram down to 10 — 5 Mo. 

all substructures that populate the galactic halo. Integrat- 
ing the conditional mass function at any redshift in different 
mass bins, the full hierarchy of substructures was automati- 
cally taken into account: subhaloes within subhaloes within 
subhaloes, and so on. The 7-ray emission process due to dark 
matter annihilation is proportional to the square of the DM 
density, and the possibility that subhaloes mi ght boost the 
expected flux was taken into consideration (see lGiocoli et all 
(|2008bl ): lpieri et all (|2008ah and references therein). 

In the present work we compute the 7-ray emission by 
considering as substructure only the progenitor haloes ac- 
creted by the main branch of our merger tree. This proce- 
dure does not take directly into account the whole 
substructure population but only the first hierar- 
chy of the population. Progressing in this way, we 
can isolate the contribution to the boost from the single 
DM structure, including its own sub-substructures subse- 
quently, following its merging history tree. Moreover, 
in this way we can model correctly the spatial satellite dis- 
tribution inside the Galaxy, since many small mass halos 
will be in fact embedded inside larger ones. 

In Fig. [2] we plot the satellite mass function, i.e. the 
mass function of progenitor haloes accreted directly by the 
main halo progenitor at any redshift. Satellite masses are 
expressed in units of the present-day mass of th e host halo 
ijvan den Bosch et al J 12009 : IGiocoli et al.1 l2008af l . This fig- 
ure shows that the mass distribution is well described by a 
single power law dA r /dln(m) oc m~ a , with slope a = —0.9, 
plus an exponential cut off at large masses. The slope was 
obtained computing the least squares fit to the data down 
to m = 10 ~ 5 Mq. This slope is slightly steeper than that 
obtained bv lGiocoli et al ] (|2008al ) using numerical iV-Body 
simulations ( « s im ~ —0.8). Such a discrepancy is proba- 




ta 



Figure 3. Satellite mass function in term of the rescaled param- 
eter f(zm, m) = S c (z m )/ yj s(m). The dotted histogram show 10 4 
realizations of the merger tree, while the solid curve is the Eq. 10. 

bly due to the fact that the present merger tree is based on 
the spherical collapse model, while numerical simulations 
are be t ter described using elli p soidal collapse ISheth et alj 
(|200ll ): ISheth fc Tormenl feoot ; IGiocoli et af] (|2007ft . 

The satellite mass function in Fig.[2]was obtained using 
the Monte Carlo code described in Sec. [2] We considered a 
halo of mass Mo at redshift z = and followed the main 
branch of its merging history tree back in time; we then: (i) 
generated a sample of progenitors at redshift 5 c ,i = <5 c ,o + 
A5 C ; (ii) identified the main (most massive) progenitor halo; 
(iii) re-run the partition code computing its progenitors at 
redshift 8 C ,2 = 5c,\ + A<5 C . We iterated these steps until the 
mass of the main progenitor halo dropped below the mass 
resolution m e . This produced the first merging history tree. 
We repeated the whole procedure, starting from the same 
mass Mo, to generate 10 4 histories, to be averaged on. 

After computing the satellite mass function (also called 
unevolved subhalo mass function) we associated to each 
halo a concentration parameter related to the quantity 
v{z m . / m) — S c (zm)/ y s(m), where z m is the satellite merg- 
ing redshift onto the host halo, and m its virial mass at z m . 

By definition, v defines the rareness of the density peak 
the halo belonged to at the epoch of merging. Higher peaks 
of density fluctuations correspond to a collapse happened 
at higher redshifts. Such a rare halo will therefore be more 
concentrated than the bulk of halos with the same mass 
whi ch formed at lat e r epoch s. 

iDiemand et al.l (|2005bh used iV-Body simulations to 
show that the present-day subhalo distribution preserves 
memory of their initial conditions. In detail, high density 
peaks are found to be more centrally concentrated and to 
move on more eccentric orbits than the overall mass distribu- 
tio n. This correlation has b een interpreted and parametrized 
bv lDiemand et al.l (|2005bh using the variable v = 8 c (z)/^/s. 
This variable (i) is related to the subhalo concentration - 
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which in turn determines the 7-ray emission, and (ii) en- 
ables one to comp ute the spatial distribu tion in the host 
halo (see Eq. (1) of IDiemand et all l|2005bl )h The details of 
the models are given in Sec. 15.21 

In Fig[3]we plot the satellite mass function in terms of 
the universal variable v. The dotted histogram shows the 
result of 10 4 Monte Carlo realizations of the partition al- 
gorithm (as in Fig[3]), while the solid curve is the Gaussian 
distribution, Eq. (jj). 

The figure shows that the satellite mass function is well 
described by a gaussian distribution for small and interme- 
diate masses, with a cutoff at u=2. This fact can be qual- 
itatively understood by recalling that the progenitor mass 
function has a Gaussian distribution when expressed in the 
rescaled variable v. For each evolutionary step dz, the satel- 
lite mass function is obtained by removing from all progen- 
itors the most massive one. Integrating over redshift, the 
total satellite mass function is effectively an integral over 
different gaussians deprived of their most massive progeni- 
tor. 



4 DRACO 

The algorithm described in the previous sections can also be 
used to study and constrain some properties of the dark mat- 
ter halos hosting the Milk- Way satellite galaxies. As an ex- 
ample let us consider Draco, a dwarf spheroidal galaxy with 
an old stellar population, dark matter dominated: Draco's 
mass to light ratio is estimated to be m/ L ~ IOOMq / Lq 
l|Prvorill992l : llrwin fc Hatzidimitrioulll995T ). This makes it a 
very interesting target to investigate the existence of dark 
matter; specifically, we study the possibility of an observa- 
tion of 7-ray due to DM annihilation inside Draco using the 
Fermi Large Area Telescope (Fermi-LAT, formerly known as 
GLAST). To this extent, we model the distribution of dark 
matter subclumps inside Draco. 



4.1 Mass and Merging Time Estimates 

Hierarchical models of galaxy formation predict that galax- 
ies are hosted by dark matter haloes. Such haloes grow along 
cosmic time through repeated merging events. In this sce- 
nario satellite galaxies reside in subhaloes accreted by the 
main halo progenitor. Satellite haloes grow themselves hi- 
erarchically, accreting mass and sub-progenitors up to the 
time when they merge onto the host halo. After they fall 
in the gravitational potential well of the host, they start 
to loose mass due to gravitational heating, tidal stripping 
and close encounters with other subhaloes, so that only a 
fraction o f their initial virial mass is still self-bound at red- 
shift zero (Ivan den Bosch et al.ll2005l ; IDiemand et al.ll2007bl . 



l200Sl ; lGiocoli et al.ll2008X 

In order to build the history tree for a subhalo hosting 
a satellite galaxy, we need to know its virial mass at merging 
time z m and use it as our starting point. From its merging 
redshift we then can go backwards in time and reconstruct 
its subhalo population, in order to model its present-day 
sub-subhalo mass function and distribution. 

Using deep wide-field multicolor CCD pho tomet r y from 
the Sloan Digital Sky Survey, assuming a iKind (|l966l ) 
spherical model of equivalent size as a reference and 
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Figure 4. Unevolved and evolved subhalo mass function of Draco 
candidates accreted by a Milky- Way halo. The solid histogram 
shows the satellite distribution, the dotted one refers to the 
present day mass fu nction evolved by us ing the average mass 
loss rate derived by iGiocoli et al . (2008a). The hatched region 
bounds the most likely values for the prese nt-day mass of Draco 
l|Odenkirchen et alJlioOll : lLokas et all2005h . 



adopting a line-o f -sight velocity dispersion of 10.7 km/s 
lArmandroff et al l (Il995h finds 3.5 ± 0.7 x 10 7 M Q within 
28 arcmin while lOdenkirchen et al.l (|200ll ). considering 8.5 
km/2, finds 2.2 ± 0.5 x 10 7 M within 40 arcmin. Differ- 
ences in the experimental mass modeling determine the 
mass range we allow for the Draco-like satellite in our 
analysis. Considering all the stars within the tidal radius, 
lOdenkirchen et all (|200ll ) determined the total luminosity 
of the Draco dwar f galaxy to be (L/L ); = 2.4 ± 0.5 x 10 5 . 
lLokas et af] |2005h also studied the distribution of dark mat- 
ter in Draco by modeling the moments of the line-of-sight 
velocit y distribution of i t s star s from the velocity dispersion 
data of I Wilkinson et all (|2004h . and obtained a best-fitting 
total mass equal to 7 x 10 7 Mq . The inferred mass-to-light 
ratio (in the V-band) was 300 Mq/Lq, almost constant with 
radius. On the other hand, we lack a direct estimate for 
the merging redshift of Draco onto the Milky Way. From 
its initial position and velocity, the accretion time of Draco 
can be derived indirectly by considering the one-to-one re- 
lation between virial radius and accret ion time implied b y 
the spherical second ary infall model (|Bertschingerl Il985l ) . 
lHavashi et all l|2003h showed that an upper limit for Draco 
merging redshift is z m < 2.8. 

We have computed the Draco virial mass from the pre- 
vious mass and redshift estimates in the following way. From 
each one of the 10 4 Monte Carlo merging trees of a Milky- 
Way sized halo (M = W 12 M /h) we noted down all satel- 
lites with mass larger than 3 x 10 7 M© accreted by the main 
branch between z — 2 and z = 2.8, totaling 8932 Draco can- 
didates. The redshift interval roughly corresponds to a time 
interval of lGyr from the upper limit of the accretion time 
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computed bv lHavashi et alj (|2003h . The lower limit on the 
Draco mass guarantees that evolved masses will lie today 
within the experimentally estimated mass range. In Fig. [4] 
the solid histogram shows the mass function of these selected 
satellite haloes. Once inside its host, each satellite will loose 
mass due to gravitational heating and tidal stripping effects. 
We modeled t hese e ffects using the results in Eq. (10) of 
iGiocoli et al.l (|2008al ). In that paper the authors measured 
the subhalo mass loss rate in a sample of high resolution 
iV-body haloes ranging from 10 115 to 10 15 M Q /h. Follow- 
ing their recipe, the ratio of z — self bound mass to the 
original virial mass of the subhalo is uniquely determined 
by the amount of time the subhalo has spent inside its host. 
The dotted histogram in Fig. [5] shows the evolved subhalo 
Draco-candidates mass function. C onsidering tha t the m ea- 
sured fractional mass loss rate by IGiocoli et al.l (|2008aT ) is 
independent of the subhalo mass, the evolved mass function 
has the same slope of the unevolved one. 

4.2 The Sub- Tree 

In order to compute the Draco satellite mass function, each 
candidate of our sample has been further evolved back in 
time along its main branch using the partition code (as it 
was done for the Milky- Way halo) considering the same mass 
and time step resolution. In this case each merger history 
tree starts at the corresponding satellite merging time. To 
increase the statistical significance of our result we run three 
realizations for each of the 8932 Draco- like satellites, totaling 
26796 Monte Carlo merger tree realizations. In Fig. [5] we 
show the unevolved sub-subhalo mass function accreted by 
the Draco sample. Comparing Fig. [5]and Fig. [5] we see that 
the two accrete d mass functions are indistinguishable: in 
fact, a s found bv lGiocoli et al.l l^OOSal h lvan den Bosch et al.l 
the accreted mass function is universal, independent 
both of final mass and observation redshift. The slope of the 
power law is again a — —0.9. 




-15 -10 -5 
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Figure 5. Satellite mass function accreted by the main branch of 
the Draco-candidate trees, until the merging on the Milky- Way 
main progenitor halo. The total number of realizations is 26796. 
Comparing the histogram with Fig. [2] we notice that the slopes 
of the dis tributions are identic a l. This is in agreemen t with what 
found bv lGiocoli et al.l ll2008ah ; Ivan den Bosch et alj J2005I") that 
is the unevolved subhalo mass function is independent both of 
final host halo mass and redshift. 

In this paper we will set $ pp = 1 and refer to $ cosmo 
as a scaled photon flux. In the following subsections we re- 
view different ways to compute <J> cosmo j n the presence of 
substructures. We will then show the computation relative 
to the merger tree technique presented in the previous sec- 
tions and will show the comparison between the different 
predictions. 



5 7- RAY FLUX FROM DARK MATTER 
ANNIHILATION IN 
SUB(SUB)STRUCTURES 

The photon flux expected from DM annihilation in the pop- 
ulation of galactic subhaloes can be modeled as 

^L(E„^AQ) = ^(E-y) x ty.An) (12) 

where ip defines the direction of observation with respect to 
the galaxy center, and AQ is the detector angular resolution. 

The particle physics dependence in Eq. [12] is given by 
the annihilation spectrum and DM properties and is embed- 
ded in the term: 

-dE^^-^^nJ'^^ 8 ^ (13) 

where m x is the DM particle mass, tr mn t; the self- 
annihilation cross-section times the relative velocity of the 
two annihilating particles, and dN^/dE-, the differential 
photon spectrum for a given final state / wi th branching ra- 
tio B f which we model after the results of ()Fornengo et al.l 
12004 ) . 



5.1 Results on the 7-ray flux from the analytical 
description of the subhalo population 

IPieri et alj (|2008al ) have considered the existence of a pop- 
ulation of substructures inside a DM halo, described by the 
mass function 

dAVdln(m) oc m~\ (14) 

and assumed that Eq. [TJ] describes the mass function of 
all substructures in our Galaxy as well as that of sub- 
substructures, at z=0. In this case all the halos that ever 
accreted onto the main one and survived are considered, 
including the higher level ones (substructures within sub- 
structures). Using this prescription it is then impossible to 
discriminate a sub-halo from a sub-sub-halo, affecting in a 
wrong way the spatial distribution of substructures inside 

the main halo. 

IPieri et all (|2008al) assumed that the substructure spa- 
tial distribution traces that of the underlying host mass from 
r V ir and down to a minimum radius, r m i n (m), within which 
disruption by Galaxy tides and stellar encounters becomes 
relevant. 
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Folding these indications together they modeled the 
number density of subhaloes per unit mass at a distance 
R from the Galactic Center (GC) as: 



Psh{m, R) = Am 



8(R- 



■H) 



(R/tf")(l + R/r? H ) 



(15) 

where rf H represents the scale radius of the parent halo 
fPH) referred to a Navarro, Frenk and White (NFW) 
l)Navarro et al.lfl997l ) profile. The effect of tidal disruption is 
accounted for by the Heaviside step function 9(r — r m in(m)). 

To determine the tidal radius, r m i„(m), they used the 
Roche criterion and compute it as the minimum distance at 
which the subhalo self-gravity at r a equals the gravity pull of 
the host halo computed at the orbital radius of the subhalo. 

To normalize Eq.[20l they imposed that 10% of the PH 
mass is distributed in subhaloes with masses in the range 
[1(T 5 - 1CT 2 ] Mph- 

Following the previous prescriptions, they found about 
53% of the MW mass (M PH = 10 12 M Q //i, r a = 21.7 kpc, 
C200 = 9.8) condensed within ~ 1.5 x 10 16 subhaloes with 
mas ses in the ra nge 

IPieri et al.l (|2008bh repeated the calculation for the 
Draco Galaxy, for which ~ 40% of the halo mass (Mnraco = 
7 x 1O 7 M , r s = OAkpc, c 20 o = 21.2) condensed into ~ 10 12 
halos with masses between [h 10~ 6 , 10 6 ]M Q /7i. 

For each substructure, they used a NFW density profile 
whose concentration parameter C200 (referred to the radius 
enclosing a d ensity equal to 200 t imes the critical density) is 
given by the lBullock et all l|200lf ) prescriptions for the sub- 
haloes at redshift zero (assuming the subhalo es are slightly 
more concentrated than the field halos). The lBullock et aD 
l|200lf ) model holds down to masses of the order of 10 5 M . 
In order to cover the whole subhalo mass range, they per- 
formed a double mass slope extrapolation down to the small- 
est masses, obtaining a concentration parameter of about 70 
for a 1O _6 M halo at redshift zero. This value is in agree- 
ment with the extrapolation a t redshift zero of what found 
in the numerical simulations bv lDiemand et all (|2005al ) , who 
isolated 1O~ 6 M halos at z=26. 

Such a model for the concentration parameter has been 
called -B- , re / model. 

The contribution of unresolved substructures to the an- 
nihilation signal along a cone of sight is given by 



<F 



dm 



"II 

J J A! 



dA 



[p s h{m, R(R Q , A, ip, 6, <£)) x P(c) x 

x <&STK c, r(A, A', V, 0', </>')) x J{x, y, z\\, B, <j>)] (16) 

where AQ is the solid angle of observation pointing in the 
direction of observation i/j and defined by the detector an- 
gular resolution 9; J(x,y,z\\,9,<fi) the Jacobian determi- 
nant; R the galactocentric distance, which, inside the cone, 
can be written as a function of the line of sight (A), the 
solid angle (8 and <j>) and the pointing angle if) through 
the relation R = VA 2 + Rq 2 - 2\R @ C, where R Q refers 
to the distance of the Sun from the Galactic Center and 
C = cos(#) cos(V') — cos{(f>) sin(<9) sin(?/>); r is the radial coor- 
dinate inside the single subhalo located at distance A from 
the observer along the line of sight defined by if) and con- 



tributing to the diffuse emission. Finally, P(c) is the lognor- 
mal distribution of the values for the concentration param- 
eters. The expression 

® C haio ° (m, c, r) = / / d<j>'d6' d\' 
J JAn ii.o.s 



p 2 (m,c,r(A,A',V,0'0')) 



A 2 



J(x,y,z\\',0'(t>') 



(17) 



describes the emission from each subhalo and p x (m,c,r) is 
the Dark Matter density profile inside the halo. 

Numerical integration of Eq. [16] gives the contribution 
to <j> cosmo from unresolved clumps in a 10 -5 sr solid angle 
along the direction tp. 

We show the results obtained in this theoretica l frame- 
work with a dashed line in Fig. [BlllPieri et al.ll2008aD for the 
MW and in Figl7l (|Pieri et alj|2008bl ) for the Draco galaxy. 



5.2 Results on the 7-ray flux from the analytical 
method including the effect of the merging 
epoch 

A furt her degree of detail was given by iGiocoli et al.l 
(20081)), who computed Eq. [TH] including the dependence of 
the subhalo spatial distribution from the initial conditions 
when the haloes accreted into the present-day Milky Way 
halo. Such a model is characterized by the u niversal variable 
u(m). According to lDiemand et aD (|2005bh the DM density 
profile of our Galaxy can be written as 



Px(r) 



(*)>(*)'] 



(18) 



with (7,/3,a) = (1.2,3,1). Subsequently, the following pa- 
rameterization is used to reflect the fact that material ac- 
creted in areas with high density fluctuations is more con- 
centrated toward the centre of the galaxy, and has a steeper 
outer slope: 



ft, = exp(i//2) 



P — > f3 v = 3 + 0.26i/ 



(19) 



Including a step function to take into account tidal disrup- 
tion, as in Eq. 1201 the number density of subhaloes per unit 
mass at a distance r from the GC, for a given v(m), becomes: 



p ah (m,r,v) = 



Am 2 6(r - r min (m)) 

G9„-7)/a ' 



)1 



(20) 



in units of M Q kpc 3 . The mass dependence in r„ reflects 
the mass depende n ce of t he virial parameter r 3 — rw/cw. 

IGiocoli et all (|2008bh normalized the number of sub- 
haloes such that 10% of the MW mass is distributed in 
subhaloes with masses in the range [10 -5 - 10" 2 ] Mmw, 
ending up with 2.4 x 10 16 subhaloes with masses between 
1O _6 M and 10 10 M Q /h, accounting for 74 % of the MW 
mass (M MW = 1.4 x 10 12 M , r a = 26 kpc). 

The contribution to $ cosmo is given by: 



<!> c 



\ip, Afi) = dm 

J M 



'"II 

J J Ai 



dOdcf, 



dA / dc 
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[p sh (m,R(R Q ,\,tp,e,<f>),v) x P{u{m)) x P(c(m))x 

x $™(™, r(\, A', ip, 9' , 4>), v, c) x J{x,y,z\X,9,<j>)] (21) 

which accounts for the influence of cosmology in the flux 
computation. P(u{m)) is the probability distribution func- 
tion for the peak rarity v(m), calculated using the ex- 
tended Press-Schechter formalism. P(c(m)) is the lognor- 
mal probability distribution for c centered on c v i r (m). While 
P(u{m)) is determined by the merging history of each sub- 
halo, P(c(m)) descr ibes the scatter in concentration for 
haloes of equal mass ()Bullock et al.ll2~00ll ) and the two prob- 
abilities may be assumed to be independent. As in the pre- 
vious analytical estimate of the subhalo population, no dis- 
tinction can be made with this method between sub-halos 
and sub-sub-halos. 

The result of this calculation for the MW is shown with 
the long-dashed line in Fig. [6] 



5.3 Results on the 7-ray flux from merger tree 
technique 

In this work we have described a merger tree approach to 
infer the subhalo population of both the MW and Draco. In 
this case, we don't need to use the P{v) and diV /dm factors, 
nor to apply any normalization. Indeed, the output of the 
merger tree gives us directly the number of objects with a 
given mass and a given v, which we call now N(fh, v). 

The total number of sub(-sub)structures found with this 
method is ~ 2.6 x 10 15 for the MW and 2.7 x 10 13 in the 
case of Draco, at the merging epoch. In the case of the 
MW, and differently from the analytical methods described 
above, this number represents only the subhaloes, while sub- 
subhaloes must be treated separately, as it has been done 
for the Draco-like subhalo. 

In order to compute the 7-ray flux we can re-write the 
subhalo distribution function as 



Psh = Ag(r,u(m))f(m) 



so that 
Ntot = 



ot = ^N{rrn)= f dV [ dm [ P{u{m))p sh 

Jn . J gal JM Jv 

which can be rearranged into the following expression: 
^N mt = I dmAf( mi ) [ P(v{rm)) f dVg (r, u(mt)) . 

m . JM Jv J gal 



Let's define 



G( mi ) = / dVg(r,v( mi )). 

J gal 

The expression of $ C osmo can be then written as: 
$ COBmo (i/>,An) = f dm [ dv [ [ d9dcj> [ d\ [ dc 

JM Jv J J AQ il.o.s Jc 



[Ag(r,v(m))f(m)P(v(m))P(c(m))$Z%°^^J]. 

With respect to Eq. [5T] we have only written p 8 h in an ex- 
plicit way and multiplied by = 1. In this way we can 



recognize and use our merger tree output ^2 N mi in the 
following way (J gal dV = J An d6d<j} J l o s dX): 

r smo (fAS!)=^ f dddtf) f dX [dc 

_. JA!1 Jl.o.s Jc 



[N(mi)g(r,u(Mi))P(c(M)) 



G(nii 



icosmo r 
-^halo ■' ■ 



(22) 



The result of this calculation is depicted with a dotted 
line in FigE] We note that the analytical methods predict 
more substructures, and therefore a higher signal, than the 
merger tree. This is mainly due to the different mass function 
slope (-2.0 versus -1.9). 

In the case of the Draco galaxy, the merger tree was 
computed at the Draco merging epoch (no sub-subhaloes 
accrete after merging), when its mass was 9.78 x 1O 9 M0. Al- 
though the sub-subhaloes evolve in redshift and loose mass 
as the parent halo does, for each sub-subhalo we consid- 
ered a NFW profile whose concentration parameter has been 
computed for the sub-subhalo mass at the merging epoch, 
motivated by the fact that the inner profile of the structures 
should not or poorly be affected by evolution. 

From the observational point of view, we are interested 
in the Draco-like galaxy as it is today, e.g. with a mass 
content reduce d by a factor ~ 10~ 2 due to tidal interactions 
after merging (jGiocoli et alJl2008al ). 

To estimate the number of sub-subhaloes in Draco to- 
day we have retained all and only the subhaloes which, at 
the epoch of merging, were inside the radius containing the 
Draco mass today, that is within r c — 2.84 kpc. r c is also 
similar t o the tidal radius ob tained using the prescription 
given bv lSpringel et all (|2008l ') to obtain the number of sub- 
subhaloes today. The total number of subhaloes is 6 x 1CF J 
smaller than the initial one. This is the upper value for the 
number of subhaloes, since we are not considering here that 
half of the subhaloes ex it the virial radius of the parent halo 
during their first orbit (|Tormen et al. I [20051 ) . In the case of 
field parent halos they would then be re-attracted inside the 
halo, but if the parent halo is a subhalo itself (like Draco) 
they would then be dispersed in the main halo (that is the 
MW). 

As for p s h, we have used the scale quantities r s and p s 
computed for the evolved Draco-like galaxy, that is for a 
7 x 10 7 M Q halo. 

The result of the computation of <& CO smo for Draco is 
shown with a dotted line in Fig[7] For comparison we also 
show the result from the analytical calculation and the same 
quantity relative to the smooth NFW halo of Draco. As in 
the case of the MW, we note that the merger tree approach 
gives a comparable yet smaller result than the analytical 

one.^ 

Pieri et al.1 (|2008af ): iGiocoli et all ( |2008bh : IPieri et al.1 



(|2008bl ) showed how the detectability of a 7-ray flux from 
such a population of sub-subhaloes with the Fermi-LAT tele- 
scope serendipitously depends on a very boosted particle 
physics contribution. The present and more accurate results 
give an even lower expected flux, thus furtherly reducing 
the hope for detection, unless an exotic boost from parti- 
cle physics enters the play. In fact, iBaltz et al.1 (|2008l ) have 
computed the map of the Fermi-LAT sensitivity to point 
sources of DM annihilations, by using the released Fermi- 
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Figure 6. Qcosmo as a function of the angle of view from the 
galactic centre, computed for the MW. 
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Figure 7. & C osmo as a function of the angle of view ip from the 
galactic centre, computed for the Draco galaxy. 

LAT response functions. Draco lies in a region of the sky 
where the 5cr detection flux above 100 MeV in 1 year of data 
taking is 0g^ r = 6 x 10~ 9 ph cm~ 2 s -1 . In order to obtain 
the afore-mentioned flux we have to multiply the particle 
physics contribution by the highest value of <E> cosmo we have 
obtained, that is to say along the line of sight pointing to- 
wards the centre of Draco. Our result is that we would need a 
boost factor of 600 (120) when using the already optimistic 
scenarios where mv = 100 (40) GeV and av = 3 x 10 -26 



In the case of Draco and for the S 20jre / model, we can 
conclude that no higher order computations are necessary 
(sub-sub-subhaloes and so on). In fact, we have computed 
the boost factor due to substructures for any mass of the 
parent halo from 10 -3 to 10 10 Mq, that is to say the ratio 
of the integral of the density squared over the whole galaxy 
including subhaloes to the same quantity computed in the 
case of a smooth halo, which gives the boost to the 7-ray 
flux due to substructures for point-like halos. We found that 
its value is pretty flat and close to 1, only slightly larger for 
larger mass haloes, which however does not appear point- 
like at the distance of Draco. This means that no boost 
is obtained when considering sub-substructures for point- 
like subhaloes. As far as massive haloes (with mass greater 
than 10 3 M Q ) are concerned, which are not point-like at the 
distance of Draco, we have to consider the fact that, as in 
the case of Draco, the boost must be distributed spatially. 
Indeed, this spatial distribution of the boost factor leads to 
an even smaller value of the expected flux along the line 
of sight toward the centre of the halo, where the flux is 
higher. This means that including sub-substructures does 
not increase the total 7-ray flux from the subhalo. 



6 DISCUSSIONS AND CONCLUSIONS 

In this paper we used the SL99 merger tree technique to 
study the mass accretion history in satellite of a Milky- Way 
sized halo. The partition code has been generalised for a 
ACDM power spectrum, considering a microsolar mass res- 
olution of mj = 1O _6 M0. The MW-main progenitor halo 
has been followed along consecutive time steps from the 
present day down to when its mass dropped below the res- 
olution limit. From the satellite haloes we have identified a 
sample of Draco-like systems and also followed them along 
their merger tree, starting from their accretion redshift into 
the MW-main halo progenitor. Both the Milky- Way and the 
Draco satellite mass function turned to be equivalent, testing 
the universality of this distribution. The partition technique 
allowed us to build up models for the present-day subhalo 
mass function for both MW and Draco - considering only 
haloes that they accreted along the main branch. 

We then computed the expected 7-ray flux from DM an- 
nihilation in such a population of substructures. We found 
that the prediction for the 7-ray flux is indeed more pes- 
simistic than the ones obtained in previous estimates which 
were not correctly taking into account the embed- 
ding of small scale sub-subhaloes within subhaloes, 
and therefore adding their contribution to the host 
halo total flux instead of to the subhalo they belong 
to. Detection with the Fermi-LAT telescope is probably out 
of the discovery range of the satellite, unless some exotic 
particle physics could boost the signal significantly. 

While we were writing this paper, the results of the 
Aquarius simulation came out (|Springel et al.ll200St ). pre- 
dicting a smaller number of subhaloes and a shallower in- 
ternal slope for the subhaloes density profiles. We will not 
repeat our analysis in their model since it would fourthly re- 
duce the expected flux, but the reader should be aware that 
the results presented in this paper, though pessimistic, are 
upper limits for the expected signal from DM annihilation 
in galactic subhaloes. 



© 2007 RAS, MNRAS OOO.nTOl 
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